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Abstract 

Tyrosine kinases (TKs) specifically catalyze the phosphorylation of tyrosine residues in proteins and play essential roles in 
many cellular processes. Although TKs mainly exist in animals, recent studies revealed that some organisms outside the 
Opisthokont clade also contain TKs. The fungi, as the sister group to animals, are thought to lack TKs. To better understand 
the origin and evolution of TKs, it is important to investigate if fungi have TK or TK-related genes. We therefore 
systematically identified possible TKs across the fungal kingdom by using the profile hidden IVlarkov Models searches and 
phylogenetic analyses. Our results confirmed that fungi lack the orthologs of animal TKs. We identified a fungi-specific 
lineage of protein kinases (FslK) that appears to be a sister group closely related to TKs. Sequence analysis revealed that 
members of the FslK clade contain all the conserved protein kinase sub-domains and thus are likely enzymatically active. 
However, they lack key amino acid residues that determine TK-specific activities, indicating that they are not true TKs. 
Phylogenetic analysis indicated that the last common ancestor of fungi may have possessed numerous members of FslK. 
The ancestral FslK genes were lost in Ascomycota and Ustilaginomycotina and Pucciniomycotlna of Basidiomycota during 
evolution. IVlost of these ancestral genes, however, were retained and expanded in Agaricomycetes. The discovery of the 
fungi-specific lineage of protein kinases closely related to TKs helps shed light on the origin and evolution of TKs and also 
has potential implications for the importance of these kinases in mushroom fungi. 
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Background 

Proteins undergo various post-translational modifications such 
as ribosylation, acetylation, thiolation, and ptiospliorylation. In 
eukaryotic organisms, reversible protein pliosphorylation acliieved 
by protein kinases (PKs) and pliosphatases plays critical roles in the 
regulation of enzyme activity and intracellular signaling. Most PKs 
catalyze ATP-dependent phosphorylation of Serine (Ser) or 
Threonine (Thr), and some of these, which are known as dual- 
specificity kinases, can also phosphorylate on tyrosine (Tyr) [1,2,3]. 
Tyrosine kinase (TK) is a distinct group that specially catalyzes the 
phosphorylation of Tyr residues in proteins. In animals, TKs play 
essential roles in cell proliferation and differentiation, immune 
responses, organ development, and other cellular processes [4,5]. 
Mutations in TK genes have been linked to various human 
diseases, such as cancer and immune diseases [6,7,8]. 

The Ser/Thr kinase catalytic domains are highly conserved and 
could be divided into 1 1 subdomains [1]. Tyrosine kinases contain 
highly conserved catalytic domains similar to those in protein Ser/ 
Thr kinases but with unique subdomain motifs. Three motifs in 
subdomain VI, VIII and XI are highly conserved in TKs but are 
not found in Ser/Thr kinase [1,9,10]. The high degree of 
conservation of the tyrosine kinase motifs could be used to 
distinguish TKs from Ser/Thr kinases. 

To date, most TKs were found in metazoan species. Previous 
studies have demonstrated that TK genes underwent duplication 



and loss during the evolution of metazoans [11]. A further study 
on dozens of eukaryotic genomes revealed that TKs appeared 
early in the common ancestor of metazoans and expanded after 
the divergence of the metazoans, especially after the split of the 
vertebrate lineage from the Ciona linage [12]. More recently, TKs 
were demonstrated to be established before the divergence of 
filastereans from the Metazoa and Choanoflagellata clades [13]. 
Many organisms outside the Opisthokont clade, such as 
Amoebozoa Acanthamoeba castellanii, Dictyostelium discoideum and 
Entamoeba histolytica [13], green alga Chlamydomonas reinhardtii [14], 
and oomycete Phytophthora infestans [15] were also found to contain 
TK or putative TK genes. 

Fungi were found to have tyrosine kinase-like kinases (TKLs) 
[16], a group of kinases that share high sequence similarity with 
TKs but function mainly as serine-threonine kinases. However, it 
is generally thought that fungi lack TKs [13,16,17]. Recently, 
possible TK genes were identified in the basidiomycete Laccaria 
bicolor using sequence searches [16] but whether these genes are 
true TKs remains to be determined. To systematically investigate 
whether TK genes occurred in fungi, in this study we searched for 
possible TKs across the fungal kingdom by using Profde hidden 
Markov models (HMMs) [17] and determined their relationships 
with TKs by phylogenetic analysis. Our results confirmed that 
fungi lack orthologs of animal TKs. However, they have a specific 
lineage of protein kinases which is most closely related to TKs. 
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Figure 1 . The distribution of FslK members in fungal species used 

of a-tubulins. The red bar indicates the number of FslK members. 
doi:1 0.1 371 /journal.pone.008981 3.g001 

Most of these genes were found in Agaricomjcetes oi Basidiomycota but 
neither in Ascomycota nor other phyla of Basidiomycota. Members of 
this lineage are predicted to have enzymatic activity but lack key 
amino acid residues that determine TK-specrfic activity. The 
evolution of members of this lineage was also addressed. 

Results and Discussion 

Identification of possible TKs in fungi 

To systematically search for the fungal possible TKs, we used 
the HMMER program [18] to search the predicted proteomes of 
84 fungi from phyla Ascomycota, Basidiomycota, Chytridiomycota, and 
^gomycota (Figure 1; Table SI) with the multi-level HMM library 
of protein kinases [17]. Only the fungal sequences designated as 
TKs (best matches) were selected and then deposited into the Pfam 
server for kinase domain confirmation. These sequences were 
further subject to preliminary phylogenetic analysis with classical 
TKs and TKLs downloaded from Kinbase (http://kinase.com/ 
kinbase/). We also included some representative fungal sequences 
classified as TKLs by our HMMER searches. In the resulting 
phylogenetic tree, fungal sequences identified as TKs were 
clustered into two distinct clades (Figure SI). One clade (fungal 
clade 2) was clustered into TKLs and was thus excluded from the 
following analysis. The other clade (fungal clade 1) was most 
closely related to animal TKs. To identify new sequences 
belonging to this clade, we built a HMM profile with the 18 
sequences of fungal clade 1 and combined it into the kinase HMM 
library to further search against fungal proteomes. The close 
relationships of newly identified sequences with TKs were also 
confirmed by the phylogenetic analysis. In total, we identified 241 
sequences from 14 fungi (Figure 1; Table S2). These sequences 
formed a distinct clade in the phylogenetic tree. We named this 
clade as fungi-specific lineage of protein kinase (FslK). 

Phylogenetic position of tiie FslK 

Because organisms beyond animals, including Amoebozoa A. 
castellanii, D. discoideum and E. histolytica, and Oomycete P. infestans 
were also found to contain TKs. We therefore performed a 
comprehensive phylogenetic analysis with selected representative 
members of FslK to determine their evolutionary relationship with 
all known TKs, by using two independent phylogenetic method- 
ologies: Maximum likelihood (ML) and Bayesian inference (BI). 
The green alga C. reinhardtii was also reported to have TKs [14]. 
We identified 16 possible TK sequences from C. reinhardtii genome 
and included them in our analysis. 

In the resulting ML and BI trees (Figure 2), the known TKs, 
including classic TKs from animals and choanoflagellates, and 
previously reported TKs from pre-opisthokont species E. histolytica, 
A. castellanii, P. infestans, formed a well-supported clade (named as 
TK clade). Three sequences of C. reinhardtii also fell into the TK 
clade. The FslK was clustered with a clade of C. reinhardtii (Cr clade 
1) and together formed a sister group to the TK clade. As we 
know, fungi are evolutionary' more close to animals than those of 
pre-opisthokont species. If fungi have orthologs of animal TKs, 
they should be clustered with them in the TK clade. In contrast, 
the position of the FslK clade suggests that orthologs of animal 
TKs were lost in fungi. 

Since the TK activity of members in Cr clade 1 is unclear, we 
do not know if the last common ancestor of both TK clade and Cr 
clade 1 has the TK activities. Therefore, whether the FslK 



in tliis study. The species tree was drawn based on the phylogenetic tree 



members have TK activity cannot be determined solely by the 
phylogenetic position. 

The members of FslK may have no TK activity 

We performed comparative analysis of TK unique motifs and 
specific residues related to TK activities in catalytic domain to 
explore whether FslK members have tyrosine catalytic activities. 
The three motifs in subdomain VI, VIII and XI are reported to be 
TK specific [1,9,10]. However, in our analysis the sequence 
pattern of the motif in subdomain X [CW(X)6RPXF] was found to 
be shared by TKs and TKLs (Figure S2) and therefore was 
excluded from our subsequent analysis. A new motif with sequence 
pattern [GXR(L/M)] in subdomain X was found to be TK 
specific (Figure S2) and was used in our analysis. Results of 
comparative analysis showed that the sequence patterns of FslK 
are obviously different from those of TKs (Figure 3). The key 
amino acid residues 'AARN' for stabilizing the relative positions of 
the substrate-binding site and the catalytic loop of TKs in motif 1 
and the first conserved proline (P) residue important for substrate 
recognition in motif 2 [10,19] were not found in FslK members. In 
addition, TKs have a glutamate (L) or methionine (M) residue in 
the fourth position of motif 3, while members of FslK have a 'P' 
residue in the equivalent site. These residues are important for the 
TK activity and are diagnostic for TKs, and the lacking of these 
residues in the members of FslK suggests that they have no TK 
activities. Members of Cr clade 1 do not contain these key 
residues, and thus may also have no TK activities. 

TKs, especially metazoan TKs, contain additional domains out 
of catalytic domains [12]. We examined additional domains in the 
FslK members. Different from the TKs, most of the FslK members 
have no additional domains; only 1 6 of the 24 1 members (Table 
S3) have additional domains which are not likely related to TK 
activities. 

These results together with the phylogenetic analysis suggest 
that the TK activity is most likely to be acquired by the ancestor of 
TK clade after it diverged from the last common ancestor of FslK 
and Cr clade 1. 

Distribution and evolution of FslK members 

Among the 241 FslK members, only 6 are from Chytridimiycota. 
All the others are irom Agaricomycotina o{ Basidiomycota. Surprisingly, 
no FslK sequences were detected in 62 ascomycetes examined. In 
Agaricomycotina, only the mushroom-forming fungi Agaricomycetes 
contain FslK members but the two TremeUomycetes, Qyptococcus 
neoformans and Qyptococcus gattii, lack any putative FslK (Figure 1). 

Phylogenetic analysis showed that FslK members were clustered 
into multiple sub-clades. Each sub-clade contains sequences from 
different species. Moreover, two distantly related sub-clades both 
contain sequences from Basidiomycetes and Chytridiomycetes 
(Figure 4). These suggest that the last common ancestor of fungi 
had possessed numerous paralogous genes from which the sub- 
clades were descended. All these ancestral genes may have been 
lost in Ascomycetes and also in Ustilaginomycotiiia and Pucciniomy- 
cotina of Basidiomycota. In contrast, Agaricomycetes retained most of 
these copies. Furthermore, lineage or species-specific gene 
duplications (gains) also have occurred in some Agaricomycetes. For 
example, the wood decaying fungus Fomitiporia mediterranea and the 
ectomycorrhizal basidiomycete L. hicolor each contains numerous 
FslK members, and many of them in each species have high 
sequence identity and are clustered together in the phylogenetic 



PLOS ONE I www.plosone.org 



3 



February 2014 | Volume 9 | Issue 2 | e89813 



Fungi Protein Kinases Related to Tyrosine Kinases 



Zal<A Dd 

MAPKKI Hs 
■MAPKKKI Hs 

CLKI Hs 



■CDK3 Hs 

■TTKI At 



y 



■MLKTAKI Hs 
■MLKTal<l Dm 
■MLKTakll Dm 

■RIPK ANKK2 Su 
■RIPK RIPK I Mm 
■RIPK RIPK3 Hs 
■RAF BRAF Hs 
■RAF RAFI Hs 
■RAF phi Dm 
■RAF lin-4S Ce 
■FLAP ksr Dm 

■TKL-TvagI 90729.m00 1 03 Tv 
■TKL-TvagI 84950.m00 1 05 Tv 
■TKL-TvagI 9409S.m00 1 06 Tv 
■TKL-Unique DDB0229955 Dd 
LRRK RocoS Dd 

Gdt Gdt9 Dd 
IRAK IRAK4 Hs 

IRAK BAKI At 
IRAK TMKl At 
IRAK pik-l Ce 
IRAK IRAKI Hs 
LISK DDB0229854 Dd 
LISK CCIG 09729 Cc 
LISK DDB0229866 Dd 
LISK LINK I Dm 
LISKTESKI Hs 



19863587 Cr 

19860786 Cr 

19868622 Cr 
1 9873084 Cr 
19865510 Cr 
19862766 Cr 
395322921 DS 



SPPG 066I3T0SP 

BDET 05232 BD 

LACBIDRAFT 335202 LB 




Out-grour 




Animals 

Chonanoflagellates 
Fungi 

Amoebozoa 

Plantae 

Oomycetes 



TKLs 



y 



Green Algae 
(Cr clade I) 



LACBIDRAFT 336006 LB 
392591399 CP 

392598148 CP 
336364380 SL 
43094 GT 

393221609 FM 
CCIG00II2TOCC 



Fungi 
(FslK) 



Amoebozoa 



Oomycetes 

Amoebozoa 

Green Algae 
(Cr clade 2) 



Src Src42A Dm 
Abl Abl Dm 
Csk CSK Hs 

UTK UTKOS Mb 

KIN 1 6 kin- 1 6 Ce 



FAK FAK Hs 

RTKF RTKF2 Mb 
RTKF RTKF3 Mb 



RTKL RTKL2 Mb 



EGFR ErbB4 Hs 
Ack ACK Hs 
I50.m00l09 Eh 
ISO.mOOIOS Eh 



o 

c 
rt 
c 
O 

U 
o8 

to 
rs 
E 
'c 

< 



kTK clade 



1 Amoebozoa 



PLOS ONE I www.plosone.org 



4 



February 2014 | Volume 9 | Issue 2 | e89813 



Fungi Protein Kinases Related to Tyrosine Kinases 



Figure 2. Phylogenetic position of the FslK. Phylogenetic trees were calculated using iVlaximum-likelihood (ML) and Bayesian inference (Bl) 
methods, respectively. Both methodologies gave similar tree topology. The tree presented here is the Bl tree. Numbers on major branches indicate 
SH-like approximate likelihood ratio test (SH-aLRT) probabilities/Bayesian posterior probabilities. Branches with Bayesian posterior probability less 
than 0.5 have been collapsed. The simple cladogram of eukaryotic groups on the top right corner was drawn according to the tree of life (http:// 
tolweb.org/tree/). Ac, Acanthamoeba castellanii; At, Arabidopsis thaliana; Ce, Caenorhabditis elegans; Cr, Chlamydomonas reinhardtii; Dd, Dictyostelium 
discoideum; Dm, Drosophila melanogaster; Eh, Entamoeba histolytica; Hs, Homo sapiens; Mb, Monosiga brevicollis; Mm, Mus musculus; Pi, Phytophthora 
infestans; Pr, Phytophthora raniorum; Ps, Phytophthora sojae; Su, Sea Urchin; Tv, Trichonionas vaginalis. For abbreviations of fungi see Table SI. 
doi:10.1371/journal.pone.0089813.g002 



tree (Figure 4), suggesting recent expansion occurred in these 
species. 

Members of FslK may have important functions in 
Agaricomycetes 

Sequence conservation of proteins is correlated witli their 
functions, and proteins with important molecular functions are 
more conserved because they are under higher selection pressure 
than those of less important ones [20,21]. Alignment of FslK 
sequences revealed that although some of them are truncated in 
the kinase catalytic domains, most of them have complete catalytic 
domains, and are highly conserved in residues required for 
catalysis, such as residues required for ATP and substrate binding 
(Figure 5). This suggests that members of FslK have catalytic 
activities of protein kinases. Carefully examined the sequences, we 
found many of the FslK members with incomplete catalytic 
domains are truncated due to sequencing gaps or wrong 
annotation, and the truncation of the others may be due to the 
sequence degradation caused by functional redundancy after 
recent gene duplication. 

We further investigated the expression of FslK genes by 
searching for corresponding EST data from NCBI and JGI 
database. Most genes examined were found to be expressed (Table 
S4), suggesting that these genes are functional in these organisms. 
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Figure 3. Comparison of sequence patterns in TK-specific 
motifs. A. Sub-domains of the protein kinase domain. Consensus 
sequences DL(R/A)A(A/R)N in subdomain VI and XP(1/V)(K/R)W(T/M)APE 
in subdomain VIII are specific to TKs. The motif in red [GXR(M/L)] was 
identified in this study. The motif CW(X)6RPXF in gray was found to be 
not specific to TKs in this study (Figure S2). B. The LOGOs show 
sequence patterns of the three motifs in each group. Red arrow heads 
indicate conserved amino acid residues that are diagnostic for TKs. 
doi:1 0.1 371/journal.pone.008981 S.gOOS 



Taken together, above evidences suggest that the FslK genes 
likely play important roles in fungi. Considering that multiple 
ancestral genes of FslK members were retained and further 
expanded in Agaricomycetes, some members of FslK most likely play 
important roles in controlling cellular development and differen- 
tiation processes specific to mushroom-forming fungi. However, 
their exact functions need to be experimentally determined. 

In summary, we systematically investigated possible TKs in 
fungi by using HMMs and phylogenetic analysis. Our results 
confirmed that fungi lack the orthologs of animal TKs. However, 
there is a specific lineage of protein kinases in fungi (FslK) which is 
most closely related to TKs. Kinases of the FslK clade lack key 
amino acid residues that determine TK-specific activities and 
therefore may not be true TKs. However, they contain conserved 
catalytic domains of protein kinases and thus are likely enzymat- 
ically active. Phylogenetic analysis revealed that the last common 
ancestor of fungi had possessed several FslK genes. The ancestral 
FslK genes may have been lost in Ascomycota and also in 
Ustilaginomycotim and Pucciniornycotina of Basidiomycota during evolu- 
tion. However, most of these ancestral genes were retained and 
further expanded in Agaricomycetes, suggesting that the FslK kinases 
possibly have important functions in controlling cellular processes 
specific to mushroom fungi. This discovery of the FslK protein 
kinases closely related to TKs helps shed light on the origin and 
evolution of TKs and also has potential implications for the 
importance of these kinases in mushroom fungi. 

Materials and Methods 

Data collection 

The predicted proteomes (listed in Table S 1 ) of fungal genomes 
used in this study were downloaded from the Fungal Genome 
Initiative (FGI) site at the Broad Institute (http:/ /www. 
broadinstitute.org), GenBank of NCBI, DOE Joint Genome 
Institute (JGI) (http://genome.jgi.doe.gov/programs/fungi/ 
index.jsf), and BluGen (http://www.blugen.org/). The predicted 
proteome of C. reinhardtii v/as downloaded from the JGI phytozome 
site (http://www.phytozome.net/) [22]. The HMM Library with 
all protein kinase domain definitions was downloaded from the 
Kinomer database (http:/ /www.compbio. dundee.ac.uk/kinomer/ 
index.html) [17]. 

Catalytic domain sequences of TKs and TKLs from Homo 
sapiens, Drosophila melanogaster, Caenorhabditis elegans, M. 
brevicollis. Trichomonas vaginalis, Tetrahymena thermophUa, 
and TKLs of Coprinopsis cinerea were obtained from the Kinbase 
(the kinase database, http:/ /kinase. com/kinbase) that collected the 
currently accepted classification of eukaryotic kinases [23]. The 
catalytic domains of putative TKs of Entamoeba histolytica were 
obtained from the Kinomer database. TKs of A. castellanii and 
three Oomycetes P. infestans, P. sojae, and P. ramorum identified 
in previous studies [13,15] were retrieved from GenBank. 

Identification of possible TKs in fungi 

The Hmmscan program in the HMMER 3.0 package [18] was 
employed to search the multi-level HMM library of protein kinases 
with each fungal proteome as queries using score of 20 as the 
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Figure 4. The maximum lil^eliKiood tree of FslK members. The phylogenetic tree was built with the kinase domain sequences using PhyML 3.1 
and was drawn using Interactive Tree Of Life Version 2.2.2 (http://itol.embl.de/#). The p-values of approximate likelihood ratios (SH-aLRT) plotted as 
circle marks on the branches (only p-values>0.5 are indicated) and circle size is proportional to the p-values. Filled red circles mark sequences used in 
the Figure 2. For abbreviations see Table SI. 
doi:10.1371/journal.pone.0089813.g004 



cutofE Only the fungal sequences designated as TKs (best matches) 
were selected and deposited into the Pfam server to confirm if they 
are kinases. These sequences were further subject to the 
preliminary phylogenetic analysis with classic TKs and TKLs to 
determine if they are closely related to TKs. 

Sequence alignment and phylogenetic analysis 

Multiple sequence alignments were performed with the PSI- 
CofFee program [24]. Alignments used for phylogenetic analysis 



were trimmed by trimAL [25] with gappyout model. Some 
sequences that are truncated due to wrong gene prediction were 
manually revised. 

Phylogenetic trees were constructed with two independent 
methods: Maximum likelihood (ML) and Bayesian inference (BI) 
methodologies. The ML trees were constructed with PhyML 3.1 
[26] using the best-fit model LG+F selected by ProtTestS [27], 
with SPRs algorithms and 16 categories of y-distributed substitu- 
tion rates. The reliability of internal branches was evaluated with 
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393215105_fM 
CC1G.04709T0_CC 
137190_GT 
336363641_iL 
CC1CJ3958T0_CC 
302692424_SCOM 



PTKccd00192 
138451_GA 
395326446.DS 
124459_HA 
393221530.FM 
66644_GT 
LACSDRAFTJ38S10JB 
SPPG,06613W^5P 
BDETJ16921_m 
U3197,GT 
393217904_fM 
392592866_CP 
CC1C_OOU2TO_CC 
30268321S.SCOM 
336369077_SL 
IACBIDRAFT_302134_LB 
2'I2223720_FP 
80700_GT 
23S584774_MP 
336372397_SL 
392589739_CP 

LAcmmrr_307628j8 

3932IS105_FM 
CClG_04709rO_CC 
137m_GT 
336363641,31 
CC1G_1395ST0_CC 
302692424_SCOM 




Figure 5. Multiple sequence alignments of representative members of FslK. The consensus sequences of Protein kinase domain (Pl<inase, 
pfam00069) and Catalytic domain of Protein Tyrosine Kinases (PTKc, cd00192) were used as references. Eleven sub-domains of FslK catalytic domains 
were shown. Conserved amino acid residues related to crystal structure and catalytic function [1,9,10] in protein kinases were indicated below. The 
default color scheme for ClustalW alignment in the Jalview program was used. 
doi:10.1371/journal.pone.0089813.g005 



SH-aLRT supports. The BI tree was constructed with MrBayes- 
3.2 [28] using mixed models of amino acid substitution with 16 
categories of y-distributed substitution rates, performing two runs 
for each of four Monte Carlo Markov Chains (MCMCs), sampling 
every 1 000th iteration over 1 . 1 x 1 0*"^ generations after a burn-in of 
101 samples. 

Examination of functional domains 

Conserved protein domains were searched in the Pfam database 
[29]at Sanger and the CDD database at NCBI (http:/ /www.ncbi. 
nlm.nih.gov/cdd). Sequence logos were generated by WebLogo 
(http://weblogo.berkeley.edu/) [30]. 

Supporting Information 

Figure SI Phylogenetic analysis of fungal sequences 
classified as TKs with those of classical TKs and TKLs. 

The phylogenetic tree was built with the kinase domain sequences 
using ML methodologies with SPRs algorithms and 16 categories 
of y-distributed substitution rates. The reliability of internal 
branches was evaluated based on SH-aLRT supports. The base 
tree was drawn using Interactive Tree Of Life Version 2.2.2 
(http://itol.embl.de/#). The p-values of approximate likelihood 



ratios (SH-aLRT) are plotted as circle marks on the branches (only 
p-values>0.5 are indicated) and circle size is proportional to the p- 
values. Sequences in fungal clade 1 and clade 2 were designated as 
TKs by multi-level HMM library of protein kinases. Other fungal 
sequences were designated as TKLs. Abbreviated species names 
are as follows: At, Arabidopsis thaliana; Ce, Caenorhabditis elegans; Cr, 
Chlamydonionas reinhardtii; Dd, Dictyostelium discoideum; Dm, Drosophila 
melanogaster; Eh, Entamoeba histolytica; Hs, Homo sapiens; Mb, Monosiga 
hrevicollis; Mm, Mus musculus; Ot, Ostreococcus tauri; Ol, Ostreococcus 
lucimarinus; Su, Sea Urchin; Tt, Tetrahymena thermophila; Tv, 
Trichomonas vaginalis; Vv, Vitis vinifera. For abbreviations of fungi 
see Table SI. 
(PDF) 

Figure S2 The identification of a new TK-specific motif 
in sub-domain X. Asterisks indicate the newly identified TK- 
specific motif [GXR(M/L)]. The previous reported motif 
[CW(X)5RPXF] shaded in gray is common in TKs and TKLs. 
(PDF) 

Table SI Information of fungal species used in this 
study. 

(XLSX) 
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Table S2 Catalytic domains of FslK kinases. 

(XLSX) 

Table S3 Additional domains of FslK kinases. 

(XLSX) 

Table S4 FslK genes detected to be expressed. 

PCLSX) 
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